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Abstract. A similarity solution for a low Mach number non-orthogonal flow impinging on a hot or 
cold plate is presented. For the constant density case, it is known that the stagnation point shifts in the 
direction of the incoming flow and that this shift increases as the angle of attack decreases. When the 
effects of density variations are included, a critical plate temperature exists; above this temperature the stag- 
nation point shifts away from the incoming stream as the angle is decreased. This flow field is believed to 
have applications to the reattachment zone of certain separated flows or to a lifting body at a high angle of 
attack. Finally, we examine the stability of this non-orthogonal flow to self-similar, three dimensional distur- 
bances. Stability characteristics of the flow are given as a function of the parameters of this study: ratio of 
the plate temperature to that of the outer potential flow and angle of attack. In particular, it is shown that 
the angle of attack can be scaled out by a suitable definition of an equivalent wavenumber and temporal 
growth rate, and the stability problem for the non-orthogonal case is identical to the stability problem for 
the orthogonal case. By use of this scaling, it can be shown that decreasing the angle of attack decreases 
the wavenumber and the magnitude of the temporal decay rate, thus making nonlinear effects important. 
For small wavenumbers it is shown that cooling the plate decreases the temporal decay of the least stable 
mode, while heating the plate has the opposite effect. For moderate to large wavenumbers, density varia- 
tions have little effect except that there exists a range of cool plate temperatures for which these distur- 
bances are extremely stable. 
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were in residence at the Institute for Computer Applications in Science and Engineering, NASA Langley Research Center, Hampton, VA 
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1. Introduction. A similarity solution for the steady, incompressible, isothermal, stagnation-point 
flow impinging on a flat plate at an angle has been studied independently by Stuart (1959), Tamada (1979), 
and Dorrepaal (1986). The problem is solved by combining an irrotational stagnation-point flow with a 
linear shear flow parallel to the plate. The stream function splits into a normal and tangential component, 
and it is discovered that the normal component satisfies the same non-linear equations which govern the 
orthogonal, stagnation-point flow. The tangential component of the stream function satisfies a linear equa- 
tion which is coupled to the solution for the normal component. The major consequence of adding a 
tangential component to the flow is the displacement of the stagnation point toward the incoming stream, 
which increases with decreasing angle. Dorrepaal also discovered that the ratio of the slope of the dividing 
streamline far from the plate to the slope of the dividing streamline at the plate is independent of the angle 
of attack. 

In the absence of a shear component, Wilson and Gladwell (1978) have shown that the incompressi- 
ble, planar stagnation flow is always stable to three dimensional self-similar disturbances. Since it is known 
from experimental observations that streamwise vortices may be formed in the region of stagnation flow 
past bluff bodies, Lyell and Huerre (1985) re-examined the planar stagnation flow problem with the aid of a 
Galerkin method. In particular, they verified the results of Wilson and Gladwell and characterized higher 
stable eigenvalue branches, showing that after the initial branch (found by Wilson and Gladwell) the other 
branches come in pairs. In addition, by a nonlinear analysis, Lyell and Huerre showed that this same flow 
is unstable for disturbances of sufficiently high amplitudes. 

Hall, Malik, and Poll (1984) extended the work of Wilson and Gladwell (1978) to include both the 
effect of crossflow in the ffeestream and the effect of suction or blowing at the wall. They showed that the 
two-dimensional flow considered by Wilson and Gladwell can be destabilized with sufficient crossflow. In 
particular, with no suction, the three-dimensional flow becomes linearly unstable at a critical Reynolds 
number (based on the crossflow velocity in the freestream) of 583.1. The effect of suction and blowing is to 
stabilize or destabilize the flow, respectively. The weakly nonlinear analysis and full numerical solution of 
the three-dimensional flow is given in Hall and Malik (1986). 

It is the purpose of this paper to explore the effects of density variations on a two-dimensional non- 
orthogonal, stagnation flow and then to study its stability to three dimensional disturbances. In Section 2, a 
similarity solution is obtained assuming small Mach number and an ideal gas. A stream function can be 
defined which is split into its normal and tangential components and is required to satisfy asymptotic condi- 
tions similar to the conditions proposed by Dorrepaal. The corresponding differential equations governing 
the temperature of the fluid and the normal and tangential components of the flow field are obtained. In 
Section 3 the formulation of the stability problem is given. Section 4 contains a presentation of our results 
of the stability problem, and conclusions are given in Section 5. 

2. The Mean Flow. It is well known that the continuity equation for an incompressible fluid can be 
exactly satisfied by introducing a stream function 

v = Vy . V = ~ Vj , (2.1) 

and furthermore, certain simple choices of the stream function satisfy the Navier-Stokes Equations. The two 
simple flows of interest here are 
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y = i3*T. y=^y 2 , 

which represent orthogonal stagnation-point flow and linear shear flow, respectively. The overbar represents 
a dimensional variable, and j3 is the dimensional rate of strain. The solution found by Stuart (1959), 
Tamada (1979), and Dorrepaal (1986) involved any linear combination of these two flows. The combination 
chosen here is 

y = 7 2 cos0 + x y sin0), (2.3) 

where 0 is related to the angle of attack 0 A by 

tan0 A =2tan0. (2.4) 

The streamlines given by the family of curves 

y = —x tan 0 ± Vjc 2 tan 2 0 + 2yji _l sec9 (2.5) 

are shown in Figure 1 for 6 = it/4. The dividing streamline is found by setting y = 0 and has slope 
m = -2 tan 0. 

Using (2.1) and (2.3) the solution to the potential flow is now given by 


U = j5 (y cos 9 + x sin 9), F=-P;ysin0, 

(2.6) 

— f F 2 v 2 ^ 


P = P~ -p^sin 2 © Y + 2 ’ 

t J 

(2.7) 

? = T 

A * OO ) 

(2.8) 


where P is the dimensional pressure found from the momentum equations, and T is the dimensional tem- 
perature assumed to be everywhere a constant. The solution for the constant density, but varying tempera- 
ture profile, has recently been given by Lyell (1990). Although the streamfunction (2.3) satisfies the 
Navier-Stokes equations, it does not satisfy the no-slip conditions at the plate y = 0. For this reason a boun- 
dary layer must exist near the plate, and in the usual fashion, is described below. The above potential flow 
field and its corresponding boundary layer is believed to have applications to the reattachment zone of cer- 
tain separated flows or to a lifting body at a high angle of attack. 

Following the analysis of Stuart (1959), Tamada (1979) and Dorrepaal (1986), solutions of the 
Navier-Stokes equations in the limit of zero Mach number are sought of the form 

U(x,r\) = xf'W) + ^g\t\), (2-9) 


p V =-/(q). 


( 2 . 10 ) 
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pOl) = 


Tin)' 


P = - ~r sin 2 9 ~=r x + -= 7 — (T|), 

2 VRe Re 


( 2 . 12 ) 


where T| is given by the Howarth-Dorodnitzyn transformation 

y 


r\ = j p(y')dy'. 
0 


(2.13) 


Here the temperature T , density p, and pressure P have been non-dimensionalized by their values at 
infinity, while lengths are referred to a characteristic length L and velocities to In addition, y and V 
have been scaled by VRe, Re being the Reynolds number of the flow. The viscosity p is assumed to be a 
function of the temperature. For small Mach number, P - 1 +yM 2 P so that P represents the pressure 
deviation from the ambient The Prandtl number Pr is taken to be a constant. The expansions (2.9)-(2.12) 
are required to match with the outer potential flow (2.6) and (2.7) as r\ ~^> and to satisfy the no slip con- 
ditions at the plate. 

Substituting (2.9)-(2.12) into the governing equations yields the following system of equations and 
appropriate boundary conditions for the normal component of the flow /, the temperature 7 \ and the 
tangential component of the flow g : 


(^rr+fr 


( rr + T sin 2 © = 0, / (0) = f (0) = 0. /' (~) = sin 0, 


(2.14) 


(fry+pr/r = o, t(o) = t p , n «o) = i. 


(2.15) 


(y- g"Y + fg" -f'g'+AT Vsin 0 cos 0 = 0, g(0) = £'(0) = 0, g ’ (°°) = ti cos 0. (2.16) 


Here T P is the nondimensional temperature of the plate, defined as the ratio of the dimensional plate tem- 
perature to that of the outer potential flow, and A = A cos 0 Vsin 0, A being the displacement thickness 
whose value will depend upon the plate temperature, Prandtl number, and the particular viscosity law under 
consideration. When T P = 1, the above equations reduce to the constant density case with A = 0.64790. 


The change of variables 


f(v\) = a F(ar|), 


(2.17a) 




(2.17b) 


, V COS0 . 

£Cn) = — T~ G(an), 

n ** 


(2.17c) 
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where a = V sin 0, reduces the above set of equations to a system which is independent of the angle. The 


equations in their simplest form are 

(-4- F"Y + FF" - (F' f + f = 0, 

T 

(- 4 - fy + PrFf' = 0 , 

T 

(- 4 - G")' + FG" - F G’ + A f = 0 , 

T 


F( 0) = F'(0) = 0. F'(o o) = 1, (2.18) 

i(0) = T p , /(oo) = l, (2.19) 

G(0) = G'(0) = 0, G'(oo) =rj. (2.20) 


Examination of the above equations reveals that the tangential component of the flow does not affect 
the normal component of the flow nor the temperature of the fluid. In fact, equations (2.18) and (2.19) are 
identical to the governing equations for an orthogonal flow. In Figure 2, F\ f, and G' are shown for the 
particular case of the linear viscosity law p. = 7\ 7> = 0.8, and unit Prandtl number (for Chapman’s linear 
law 1 1 = C f y the constant C can be scaled out by rescaling q, F, G, and A appropriately, but must be 
borne in mind when transforming the variables back to their dimensional forms). 

A quantity of importance is the viscous displacement of the stagnation point. Note that as T| 0, 

F=y V--yTl 3 + •••, G = y T » 2 "l L||3+ (2.21) 

where C and D are chosen numerically so as to satisfy the boundary conditions at infinity. For the constant 
density case, 7 > = 1 , C = 1.232588, and D = 1.406544. C\ and D x are known constants found by substi- 
tuting the expansions into the original equations. The viscous displacement of the stagnation point can be 
found by writing (2.21) in terms of / and g, 

/ = — q 2 (sin0) 3/2 - -^-q 3 sin 2 0 + *■*, g=-“Tl 2 cos0 — -^-q 3 cos0 Vsin0 + ■ • • . (2.22) 
2 6 2 3 


This leads to the following expansions of the velocity field near the plate 


U = x (C r| — ^ r\ 2 Vsin0 )( sin 0 ) 3/2 + -^=(D T|cos0 - Z^t^cos© Vsin 0), 


(2.23) 


V = - -y T P ( sin 0 ) 3/ V • (2.24) 

The dividing streamline near the plate is found by setting (2.23) to zero; then, taking the limit r\ — > 0 pro- 
vides the following expression for the viscous displacement of the stagnation point 

X sp = - -j= y cos 0 ( sin e )~ m . (2.25) 

For the constant density case, the ratio D / C = 1.141131. The stagnation point shifts in the direction of the 
incoming flow, and this shift increases as the angle decreases from n/2. In Figure 3 the ratio D /C versus 
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Tp is shown using the linear viscosity law and unit Prandtl number. Note that a critical plate temperature 
exists around 7.4 for which the ratio vanishes. For plate temperatures above this value the stagnation point 
shifts away from the incoming flow as the angle is decreased. 

3. Formulation of the Stability Problem. The problem considered here is the stability of the non- 
orthogonal stagnation flow described in the previous section. Rescaling the transverse coordinate z and velo- 
city W in the same manner as y and V, respectively, and keeping only the leading order terms in a large 
Reynolds number expansion, it is easily shown that the non-dimensional equations which govern the stabil- 


ity of this flow are given by 

Pi + (pU) x + (pV\ + (p W\ = 0, (3.1a) 

p [U t + UU X + VU y + WU Z ] + P X = at U y ) y + ftt U z ) zy (3.1b) 

p[V t + UV X + W y + WV Z ] + R P y = Oi U y ) x + (2pV y ) y + + H%)], - (3.1c) 

p [W t +UW x +VW y + WW Z ] + RP t = (1 \iU t ) x + (2pW 2 ) z + [p(V z 4- lVy)] y - K tJ (3. Id) 

p [T t + UT X + VT y + WT Z ] = Pr~ ] (pj, ) y + Pr~' (pT, \ , (3. le) 

pT = 1. fi = 7\ (3.1f,g) 


2 

where A' = — |i(i/ z + V y + tt^). The linear viscosity law (3.1g) has been used, and the Mach zero limit has 

been taken. Note that the pressure can be eliminated from (3.1c) and (3. Id) in favor of the vorticity com- 
ponent, Q = Wy - V 2 . 

Guided by the fact that the incoming angle can be scaled out of the mean flow equations, solutions 
are sought in the following form 

U(x,t\j ,0 = a 2 x F'(a t)) + a" 1 Re" 1,2 cos6 G'(a tj) + e lx a 2 Ul ( a rtf + R e m u 2 {a tj)] EXP , (3.2a) 


V(x,T\j,t) = -a f (a r\) F(a r\) + ea f (a t\) v t (a r\) EXP , (3.2b) 

W(x ,r\,z j) = za w\(cn\)EXP , (3.2c) 

T(xs],z,t) = f (a Tl) (1 + zT,{a r\)EXP), (3.2d) 

p (x ,n,z ,/ ) = f~\a TO (1 - e r, (a T)) EXP ), (3.2e) 

SM.x,f\j,t) = ea 2 f~\a n)Qj(a tj) EXP, (3.20 



where 
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EXP = exp{ i a az + a 2 tot ), (3.2g) 

and £<1. Note that the size of the pressure disturbance must be 0(1/Re) so as to balance equations (3.1c) 
and (3. Id) and hence does not effect (3.1b). Here, a is the spanwise wavenumber, taken to be real, and co is 
the real temporal growth rate. The assumption of co being real is consistent with previous investigations (eg., 
Wilson and Gladwell, 1978) and was verified numerically by allowing co to be complex. In every case, con- 
vergence resulted in co being real. Hall, Malik and Poll (1984) obtained complex values, but only with 
crossflow. Presumably, the effect of crossflow is to destroy the two-dimensional symmetry of the problem, 
thus allowing for complex eigenvalues. The disturbance is considered amplified, neutral, or damped 
depending on whether co is positive, zero, or negative, respectively. These stability results are limited to the 
special class of self-similar disturbances considered above. Recently, Brattkus and Davis (1990) considered a 
wider class of disturbances for the stability of the incompressible Hiemenz flow and concluded that these 
modes were more stable as compared to the self-similar disturbances. Therefore, the use of self-similar dis- 
turbances is justified for the constant density case, and there is yet no reason to believe that this will not 
also be true for the variable density case, although this still needs to be verified. 

Substituting the expansions (3.2) into (3.1) yields the linearized equations for the perturbation ampli- 
tudes 

u x + v \ + i aw { = coTj - F T x \ (3.3a) 

u +F - (a 2 f 2 + co + 2 F' ) k ] =F" - [2(F') 2 - 2F F" - f]T l - F" T x \ (3.3b) 

£}[ = - i af 2 vi, (3.3c) 


i2,"+ (F -//)£!,'- [a 2 f 2 + (t)-F +F//JQ, = <&, 


(3.3d) 


1 

Pr 


1 


+ (F 


2 t'2 


^ j_r \ rp , r Ct 1 

* H)T ' ~ [ ~p7 


+ 0) + 2 FH]T 1 =Hv l , 


(3.3e) 


where 


<!> = //[ 


■]- 


i a F Z F (1 + 2Fr) V] + i a f z u, - tow, I - 2 i af 2 F' T\ 


iaf 2 [ 


F" - F F' + H [2fO + 4F 7 - F 2 (l + 2 Pr)] 


] 


(3.4) 



(3.5) 


Note that the above system is independent of the angle 0 of the mean flow, and hence the stability problem 
for the orthogonal case is identical to the stability problem for the non-orthogonal case. The above system 
is eighth-order for the five unknowns u u v u T u and In addition, there is an equation for u 2 which 
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is non-homogeneous and involves G . Once the system (3.3) has been solved, u 2 can be found directly. 

The appropriate boundary conditions at the plate are 

«i = v 1 = k> 1 = 7’ 1 = 0. (3.6) 

Before the boundary conditions at infinity are considered, it is convenient to first rewrite (3.3) as a real sys- 
tem. To this end redefine u l =Imag(u l ), v, = Imag(y x ), w x =Real(w{), 7, = Imag{T v ), and 
Oj = Real ( Qj), where Real and I mag mean the real and imaginary parts of the complex function are to be 
taken. Then (3.3) becomes 

«i + v t ' + awj = coT, -F 7Y, (3.7a) 

«i" + F M,' - (a 2 f 2 + (0 + IF' )mj = F" vj - [2(F') Z - 2F F" - f]7, - F” 7,\ (3.7b) 


£2i = tV]' + a7 2 v lt 


(3.7c) 


0," + (F -H)G i'- [a 2 7 2 + co - F' + F H ]Q, = d>. 


(3.7d) 



ft 

1 


a 2 r 2 


+ <F + Tr mr<-i Pr 


+ (0 + 2F H] T x = // Vi, 


(3.7e) 


where 


d> = -// 


+ a 7 


af 2 / 7 (l+2/ 5 r)vi + a'f 2 Mi + C0 H'ij + 2 a f 2 F' 7/ 
F" - F F' + H [2 (a + 4 F" - F 2 ( 1 + 27r)] ] 7,. 


The asymptotic behaviour of the perturbation amplitudes at infinity are given by 

7" i = Xo ( r| — y a2/Pr * “1 + e.s.t . , 


(3.8) 


(3.9a) 


m, =X,(T1-/1 ) ( “ 2 + “ + 2) (l + 0(1)) + e .s.t., (3.9b) 

where e.s.t . stands for exponentially small terms as T| —» v lf and wi have similar algebraic terms in 

their asymptotic expansions; in addition, v ^ and W\ have an exponentially growing term e ai \ Wilson and 
Gladwell (1978) have shown that the correct boundary conditions for the stability problem is that all pertur- 
bation quantities must decay exponentially at infinity. Thus, for a fixed wavenumber a, the eigenvalue prob- 
lem is to solve the eighth-order system (3.7) subject to the boundary conditions (3.6) while determining a 
value of co so that all algebraic and exponentially growing terms are zero. 

Note that the above system is identical to that considered by Wilson and Gladwell (1978) in the con- 
stant density limit. They solved this problem by use of invariant embedding and Riccati transformations. 
Lyell and Huerre (1985) solved the same problem, but using a Galerkin method. It should be emphasized, 
however, that their main interest was the nonlinear stability. Here, the above eigenvalue problem is solved 
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by embedded shooting methods using a Runge-Kutta scheme with variable step size, which proved to be 
very fast and very inexpensive. For a given value of a, an initial © is chosen and the system (3.7) is 
integrated subject to the conditions at the plate (3.6) and 

ft, = C,. = C 2 , T\ = C 3 , Q'i = C 4 . (3.10) 

Then, an iteration on © and the C ; ’s by a Secant method is performed until X®, X 1( and similar constants for 
the other quantities were less than 10 -6 . All calculations were done in 64-bit precision, taking approximately 
one second of CPU time on the CRAY-2 for each eigenvalue. 

4. Results. In all calculations the Prandtl number is taken to be one. 

As discussed in the Introduction, for the constant density 7> = 1 case, Wilson and Gladwell (1978) 
have shown that the planar stagnation flow is always stable (© < 0) to three dimensional self-similar distur- 
bances. In their paper a table is given for selected values of the wavenumber a with corresponding values of 
the growth rate ©. In particular, the authors showed that a maximum exists in © for an a of about 0.298. 
In addition, they suggested that there were other more stable branches than were shown, but this was not 
pursued in detail. Lyell and Huerre (1985) re-examined the planar stagnation problem with the aid of a 
Galerkin method. They verified the results of Wilson and Gladwell to within 10% using a one term Galer- 
kin expansion or to within 5% using a two term expansion and displayed two more eigenvalue branches. 
The authors showed that after the initial branch, found by Wilson and Gladwell, other stable branches come 
in pairs. To verify the present numerical scheme, the first three eigenvalue branches are plotted in Figure 4. 
Results for the first branch agree with those of Wilson and Gladwell to the figures they gave. The other two 
branches, which occur as a pair, agree with the general behaviour given by Lyell and Huerre. However, 
exact agreement is not possible since a Galerkin scheme will degenerate for higher modes unless sufficient 
number of expansion terms are kept. The calculations presented here are believed to be correct to six digits. 

The purpose of this section is to extend the above work so as to include the effect of variable density. 
Results for the least stable branch are plotted in Figure 5 which shows the decay rate - © versus the plate 
temperature 7> for various selected values of a. In Figure 5a, small values of wavenumber are considered 
(a < 0.4). Note that for a = 0.001, the least stable branch ceases at some critical plate temperature slightly 
greater than one. As a increases, so does this critical plate temperature. Most importantly, note that the 
effect of cooling the plate is to reduce the value of the least stable eigenvalue, while heating the plate is 
seen to have the opposite effect; i.e., increasing the stability of these modes. If the plate temperature is 
increased much above the value of one, then each of these small wavenumber modes vanish in favor of a 
more stable branch. The plate temperature at which these small wavenumber modes vanish increases for 
increasing wavenumbers. In Figure 5b, moderate to large wavenumbers are considered (a > 0.5). For 
a = 0.5 the least stable branch ceases at some critical plate temperature less than one, about 0.3. This criti- 
cal value increases as a increases. Therefore, it is seen that the variation in density has little effect on 
moderate to large wavenumber disturbances except that there exists a range of cool plate temperatures for 
which these disturbances are extremely stable. 

It is now clear from Figures 5a and 5b that eigenvalues on the least stable branch can cease to exist 
for a particular combination of a and 7> . In order to investigate this behaviour, the dispersion relation is 
shown in Figure 6 as a sequence of graphs for a = 0.9 and decreasing values of the plate temperature. 
These graphs show the value of X, versus © after eliminating all other algebraic and exponentially growing 
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terms. Figure 6a shows the dispersion relation for the case 7> = 1. First note that for co > 1.5 the curve is 
increasing and approaches zero only in the limit to «>. Thus, X x never vanishes, and the flow is linearly 
stable. For to < 1.5, the first zero of is at co = -1.853386, which agrees with Wilson and Gladwell (1978). 
For values of co less than this, the zeros of X l occur in pairs, as was shown by Lyell and Huerre (1985). 
Once density variations are taken into account, no matter how small the change in plate temperature, the 
structure shown in Figure 6a is completely altered. Figure 6b is a graph for T P = 0.9. Notice that the verti- 
cal asymptotes have changed and additional ones have been added from Figure 6a to 6b. Also the ordering 
scheme of Lyell and Huerre for the constant density case is no longer valid. Figure 6c shows that some of 
the vertical asymptotes disappear while others merge, thereby forming a continuous curve in over a broad 
range of co. This pattern continues through Figures 6d and 6e. The important thing to note is that between 
Figure 6d and 6e, the least stable mode has disappeared. 

The eigenfunctions for the first three modes of Figure 6c are shown in Figures 7a,b,c, respectively. 
The first two modes form a pair, and hence their eigenfunction structures are similar. The third mode is iso- 
lated, and its eigenfunction structure is different from those of the first two. Also note that each eigenfunc- 
tion decays exponentially at infinity. 

The effect of the angle 0 of the mean flow on the linear stability can now be assessed from (3.2g). 
For any given value of 7> , as the angle is decreased so is the wavenumber and the magnitude of the tem- 
poral decay rate, thus making nonlinear effects more important. 

4. Conclusions. It has been found that variations in density can have a significant effect on both the 
structure of a low Mach number, non-orthogonal stagnation flow and its linear stability to three dimensional 
self-similar disturbances. As in the constant density case, the tangential component of the flow can be 
decoupled from the normal component. It is shown that the viscous displacement of the stagnation point 
due to the non-orthogonal flow is 0(Re ~ v 2 ), which depends on the plate temperature. In particular, a critical 
plate temperature exists above which the stagnation point shifts away from the incoming flow instead of 
toward the flow as in the constant density case. 

For the linear stability, a proper scaling can be defined so as to eliminate the angle of the mean flow 
when considering the stability characteristics. Therefore, the stability problem for the orthogonal problem is 
identical to the stability problem for the non-orthogonal problem. We show that, for small wavenumbers, 
the temporal decay rate of the least stable mode is decreased by cooling the plate. Finally, we show that 
decreasing the angle of attack decreases the magnitude of the temporal decay rate, thus making nonlinear 
effects more important. 

We note that the above results are restricted by the fact that the similarity solution of the mean flow is 
not appropriate for all values of x . We are currently extending the above analysis to a fully compressible 
subsonic flow by means of the IHingworth-Stewartson transformation. 
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Figure 1. Plot of the outer potential streamlines for 0 = n/4. 


- 12 - 



Figure 2a. Plot of F' versus T] for T P - 0.8 and unit Prandtl number. 
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Figure 2b. Plot of T versus t\ for Tp =0.8 and unit Prandtl number. 
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Figure 4. Plot of - to versus a for 7> = 1.0 for the first three eigenvalue branches. 



Figure 5a. Plot of - to 








for a = 0.9 and Tp = 0.9. 



a = 0.9 and 7> = 0.8. 













